Reionization history and CMB parameter estimation 
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We study how uncertainty in the reionization history of the universe affects estimates of other 
cosmological parameters from the Cosmic Microwave Background. We analyze WMAP7 data and 
synthetic Planck-quality data generated using a realistic scenario for the reionization history of the 
universe obtained from high-resolution numerical simulation. We perform parameter estimation 
using a simple sudden reionization approximation, and using the Principal Component Analysis 
(PCA) technique proposed by Mortonson and Hu. We reach two main conclusions: (1) Adopting a 
simple sudden reionization model does not introduce measurable bias into values for other parame- 
ters, indicating that detailed modeling of reionization is not necessary for the purpose of parameter 
estimation from future CMB data sets such as Planck. (2) PCA analysis does not allow accurate 
reconstruction of the actual reionization history of the universe in a realistic case. 



I. INTRODUCTION 

Temperature and polarization power spectra of cosmic 
microwave background (CMB) currently offer the most 
precise test of cosmological parameters. However the 
CMB is subject to cosmological foregrounds. If there 
is a phase of the history of the universe that we do not 
understand completely, it may bias our cosmological in- 
terpretation of the CMB. One such a phase is the epoch 
of reionization (EoR). Given the enhanced precision of 
upcoming Planck data, understanding and quantifying 
this bias is of significant importance. List of previous 
studies on this topic includes but not limited to 1 (i . 

Reionization of neutral hydrogen is the last major ther- 
modynamic phase transition in the history of the uni- 
verse. This era begins by the emergence of the first lu- 
minous sources such as early star-forming galaxies and 
rare quasars and it is believed to be completed a hun- 
dred million years later. As the luminous sources start 
to produce energetic photons, they reionize the neutral 
hydrogen around them. Free electrons generated during 
reionization will Thompson scatter the CMB photons in 
the period after recombination. This scattering results in 
small suppression in CMB temperature anisotropy on all 
scales and also will generate an extra polarization signal 
on large angular scales. A commonly used parameter to 
characterize these two effects is the total scattering op- 
tical depth r. The current constraint on the value of r 
from WMAP7 [7] derived from the measurement of the 
amplitude of E-mode of polarization is r = 0.087± 0.017. 

The above constraint on r is obtained assuming that 
reionization is a sharp transition corresponding to red- 



shift of z r 



10.5 ± 1.2. However a sudden reioniza- 



tion requires an infinite photon production rate. Since 
the galaxy mass function evolves smoothly, reionization 
should be a gradual process which takes of order of Hub- 
ble time. Therefore this unphysical assumption about 
the reionization history can introduce a bias in the con- 
straints on r and other cosmological parameters degen- 



erate with it. Another possible source of bias arises if 
temperature and polarization spectra of CMB are sensi- 
tive to details of reionization history beyond the optical 
depth. In this case, even if modeling reionization as a 
sudden process gives the correct estimate of r, inaccu- 
rate modeling of redshift evolution of number density of 
free electrons in intergalactic medium, x e , would result 
in an additional bias. 

In recent years there have been a significant improve- 
ments in numerical simulations and semi-analytical mod- 
els of reionization which provide us with a more realis- 
tic picture of this era. In the first part of our analysis 
we adapt the result of one of these simulations to in- 
vestigate the interplay between reionization history and 
CMB parameter estimation. Although the combination 
of simulation and semi-analytical models provide us with 
a more physical picture of EoR, the reionization process 
is complex, and there is currently a lack of direct obser- 
vational data on the EoR. Therefore, an accurate model- 
independent approach to reconstruction of the reioniza- 
tion history is also of great importance. The second 
part of our analysis is focused on the study of accuracy 
of the commonly used model-independent approach (Hu 
& Holder, Mortonson & Hu) for Planck-level sensitivity 
data. 

The outline of the paper is as following: in Section |TT] 
we discuss the more realistic reionization model that we 
have adapted in our analysis, in Section III we discuss 



the analysis methods while in Section |IV| we present our 
results. In Section W\ we draw our conclusions. 



II. CONSTRAINTS ON REIONIZATION 
HISTORY FROM SIMULATIONS 

We adapt the reionization history from Volonteri & 
Gnedin [8] as a physically motivated reionization model. 
The ionized fraction function of redshift for this 

model is shown in Fig. [T] This result is based on theoret- 
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FIG. 1: Reionization history from Volenteri & Gnedin [8] 



ical modeling of sources of ionizing photons and count- 
ing the number of ionizing photons per atom, while ac- 
counting the photon loss due to recombination in an ap- 
proximate manner. The photon budget in this model 
is supplied by two main sources: stellar population and 
quasars. When considering the contribution of quasars 
to the total reionizing photons, the "secondary reioniza- 
tion" by energetic photons from quasars is also included. 
Other more exotic sources such as energetic photons and 
electrons from dark matter annihilation (see [5] and ref- 
erences therein) are neglected. 

The contribution from quasars to the total ionizing 
photons is estimated by modeling the evolution of mas- 
sive black holes within dark matter halos. The contri- 
bution from the star-forming galaxies to the total ion- 
izing photons is calculated by extrapolating of the ob- 
served UV luminosity functions of high redshift galaxies 
(Bowens, el oi.[TD], [H]) and using an estimate of rela- 
tive escape fraction of the ionizing radiation from Gnedin 
(2007) PJ]. We refer the reader to Volenteri & Gnedin 
[8] for the details of this calculation; here we only list the 
main ingredients of the model for the sake of complete- 
ness. 

The total number of ionizing photons per atom is given 
by: 



N. 



j/a 



N, 



' ■-■* + + fsi(x) —j^ . (1) 



The first term in the above equation accounts for the 
contributions from stars. The second term is the con- 
tribution from ionizing photons from quasars with the 
mean energy of E 1 while the last term includes the sec- 
ondary ionizations from energetic photons. Quantity 
n a = (1 — 0.75Y p )n b ps 2.0 x 10~ 7 cm -3 is the comov- 
ing number density of hydrogen or helium atoms, 14.4 
cV is the mean ionization energy per atom, and U y is 
the total comoving energy density of ionizing radiation 
emitted by the population of quasars. The mean photon 
energy is a parameter of the model and is assumed 



to be 300eV. The exact value of E 1 depends on the spe- 
cific shape of the average quasar spectrum. However, 
since most of the contribution from quasars comes from 
secondary ionizations, which are independent of E^, the 
result is only mildly affected by the actual value of that 
parameter (see [8] for more detailed discussion). For the 
uniformly ionized gas, /si is the fraction of radiation en- 
ergy density going into secondary ionizations, which is 
calculated in [T^j and fitted by a simple formula in [2] , 

f 8I (x) « 0.35(1-^ 4 ) 18 -1.77 (^y\ - 2 (l-x°-r. 

Uniform ionization is a simplistic assumption. Therefore 
the quantity x is an effective value x c g such that 



(3) 



and the average is mass- weighted. A complete calcula- 
tion of x e fi requires precise modeling of transfer of radi- 
ation throughout the universe. In the absence of such 
modeling, Volontcri, et al. assumed the following ansatz: 



£ c ff = min(l, f n ^ x N j/a ), 



(4) 



where / n _> x is a free parameter that accounts for loss 
of ionizing photons due to recombinations of previously 
ionized atoms. In reality it may be a function of time; 
it is unlikely that / n _*. x is much smaller than one, or the 
universe would not be fully ionized by z m 6. For the 
reionization history shown in Fig. [T] f n -yx — 0.75 (our 
fiducial value). 



III. ANALYSIS 

We investigate the dependence of constraints on cos- 
mological parameters obtained from CMB on the as- 
sumed reionization history. We consider three reioniza- 
tion histories: instantaneous reionization, a physically 
motivated reionization history described in Sec. [^(Here- 
after we refer to this reionization history as VG model) 
and a general reionization history parameterized in terms 
of its principal components (PCs) with respect to E-mode 
polarization (Hu & Holder, Mortonson & Hu). We per- 
form three sets of Markov Chain Monte Carlo (MCMC) 
analyses using the publicly available CosmoMC code [TS] 
as modified by Mortonson & Hu [TB] and compare the 
constraints on cosmological parameters obtained for the 
above reionization profiles. 

To consider a general reionization history, the ionizaed 
fraction, x e is considered as a free function of redshift 
and is decomposed into its principal components (PCs) 
with respect to E-mode polarization of the CMB: 



x e (z) 



-fid 



]rriiSi{z). 



(5) 



Si are the eigenvalues of the fisher matrix that describes 
the dependance of Cf E on x e (z) and n%i are the ampli- 
tudes of the principal components. x[' ld is the fiducial 
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FIG. 2: Marginalized constraints in the <jg,-n 3 plane for (a) WMAP7 data and (b) simulated Planck data sets. Shaded blue 
contours correspond to PC reionization while dotted contours are obtained assuming VG reionization history. In (a) the shaded 
red contours correspond to sudden reionization while in (b) the constraints for sudden reionization are shown with solid red 
contours. 



model at which the Fisher matrix is calculated. The 
principal components are defined over a limited range 
in redshift, 6 < z < 30, with x e = at redshifts 
z > 30 and x e = 1 at redshifts z < 6. In our MCMC 
analysis we use the five lowest-variance principal compo- 
nents constructed around a constant reionization history 
of x{ id = 0.15 and vary the amplitudes of these com- 
ponents. We impose the same priors on m as those in 
Mortonson fc Hu [TO]. 

The data sets used are temperature and polarization 
data from 7-year WMAP (WMAP7) 17J and simulated 
Planck-precision data described below. We consider a 
flat ACDM cosmology with the base cosmological pa- 
rameters chosen to be the baryon and cold dark matter 
densities, f^/i 2 and fl c h 2 , the ratio of the sound horizon 
to the angular diameter distance at decoupling, S and 
the spectral index of scalar perturbations, n s . Follow- 
ing Mortonson and Hu, instead of varying A s , we take 
a combination of 10 10 A s e~ 2T as a free parameter. This 
removes the known degeneracy between the amplitude of 
scalar perturbations and the optical depth. The optical 
depth, t, is then not a free parameter but is a derived 
parameter. For the general reionization parameterized 
by principal components, the five lowest variance ampli- 
tudes of PCs are varied. For WMAP7 data, the priors 
on cosmological parameters are set to the best fit values 
from the WMAP team. To have self-consistent compari- 
son, for Planck data, the priors on the base cosmological 
parameters are set to best-fit values of our WMAP7 run 
with VG reionization history. Chain convergence was de- 
termined with the Gellman-Rubin statistic, R — 1 < 0.1. 

Planck data is simulated for a fiducial model with the 
best fit values from our WMAP7 run with VG reioniza- 



tion history with an optical depth of r = 0.0408. Since 
for the simulated data we know the true reionization his- 
tory for which the data is produced, we can study with 
more clarity the possible bias in parameter estimation 
that is introduced by incorrect assumption about the 
reionization history. The data is simulated in three chan- 
nels with frequencies (100 GHz, 143 GHz, 217 GHz) and 
noise levels per Gaussian beam (cj ix ) 2 = (46.25 /xK 2 , 
36 fiK 2 , 17.6 fiK 2 ) (with ct£ x = \/2crJ ix for polarization 
spectra). The Full Width at Half Maximum (FWHM) 
of the three channels are #f w hm = (9. 5', 7.1', 5.0') and we 
assume a sky coverage fraction of 0.65 [18]. 

IV. RESULTS 

We present the results of our MCMC analysis as 2- 
dimensional contour plots on base cosmological parame- 
ters. The constraints on derived parameter erg and on 
cosmological parameters for WMAP7 data are shown in 
Figs. [2] (a) and [3] respectively. The corresponding con- 
straints for simulated Planck data are shown in Figs. [2] 
(b) and|4] The dotted black lines correspond to VG reion- 
ization, while the filled blue contours are obtained as- 
suming the PC reionization. In WMAP7 plots, the filled 
red contours correspond to sudden reionization while in 
Planck plots, the corresponding constraints are shown 
with solid red contours. 

For WMAP7 data, assuming VG reionization history 
results in a smaller value of r compared to the sudden 
and PC reionizations. Due to the degeneracy between 
r and A s , a lower value of r results in smaller values of 
A s . Our results indicate that for WMAP7, constraints 
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FIG. 3: Marginalized constraints on cosmological parameters for WMAP7 data set. Shaded blue contours correspond to PC 
reionization while shaded red contours correspond to sudden reionization and dotted contours are obtained assuming the VG 
reionization history. 




FIG. 4: Marginalized constraints on cosmological parameters for Planck precision CMB data set. Shaded blue contours 
correspond to PC reionization while solid red contours correspond to sudden reionization and dotted contours are obtained 
assuming the VG reionization history. 



5 



Reionization Histories 



— Sudden 




5 10 15 20 

z 

FIG. 5: Reconstructed reionization history from the first five 
principal components using simulated Planck data shown in 
green while sudden and VG reionization histories are shown 
in red and blue correspondingly. 

on base cosmological parameters and the derived param- 
eter a% are mildly affected by the assumption about the 
reionization history (the shift is less than la). Assuming 
a general reionization history and truncating the expan- 
sion at the fifth principal components of x e , degrades the 
constraints on cosmological parameters. Nevertheless the 
constraints are consistent with those obtained assuming 
sudden reionization or the VG reionization. 

One would expect that the enhanced precision of 
Planck data would distinguish between different reion- 
ization models since they would have an impact on the 
constraints on cosmological parameters. However our re- 
sults indicate that for reionization histories that are not 
very different than instantaneous reionization, even at 
the sensitivity level of upcoming Planck data, parameter 
estimation from the CMB is only affected by the value 
of optical depth and not the details of the reionization 
history. Figures [2] (b) and [4] show the constraints on cos- 
mological parameters for three reionization histories from 
simulated Planck-quality data. The black dotted con- 
tours show the constraints assuming the true reionization 
history (VG reionization) for which the data is simulated. 
These results indicate that even at the level of Planck 
data, constraints on cosmological parameters are not af- 
fected by the assumed reionization history. Comparison 
between the filled blue contours and solid red contours 
indicates that even for Planck data, the constraints ob- 
tained assuming sudden reionization are as good as those 
obtained assuming general reionization. 

Moreover, our results indicate that principal compo- 
nent analysis does not allow for accurate reconstruction 
of reionization history of the universe from CMB data. 
Fig. [5] shows the reconstructed reionization history from 
Planck data using PCs in addition to the corresponding 
sudden and VG reionization. The reconstructed reion- 



ization using PCs exhibits oscillatory behavior at z > 10 
while the true reionization (VG) is a a monotonically de- 
creasing function of redshift. 



V. CONCLUSIONS 

The epoch of reionization is the last major phase tran- 
sition in the history of the universe, which sets a screen 
between us and CMB. Due to the complexity of the reion- 
ization process and lack of direct observational data, in 
analyzing CMB data, it is commonly assumed that reion- 
ization is a sharp transition occurring at a given redshift. 
Given the enhanced precision of upcoming Planck data, 
understanding the impact of this unphysical and simplis- 
tic assumption is of great relevance. 

In this paper we investigated the question of whether 
the assumption of sudden reionization introduces any 
bias into constraints on other cosmological parameters. 
We conducted three sets of MCMC analysis assuming 
three different reionization histories: sudden reionization, 
a physically motivated reionization history (Sec. [Tl]), and 
a general reionization parameterized by principal compo- 
nents (Hu & Holder, Mortonson & Hu) we compared the 
constraint on cosmological parameters obtained for these 
three reionization histories. Our results can be summa- 
rized as follows: 

• For reionization histories which are not drastically 
different from instantaneous reionization, such as 
the physically motivated model from Gnedin, et 
al. simulations (Figfl]) , constraints on cosmological 
parameters are not sensitive to details of reioniza- 
tion history (Figs. [3]-|4]). Therefore lack of precise 
knowledge of true reionization history and assum- 
ing sudden reionization model does not bias param- 
eter estimation from CMB. 

• Principal component parametrization of reioniza- 
tion history does not allow for accurate reconstruc- 
tion of ionized fraction. 

This result is fully consistent with previous results in the 
literature, in particular Columbo and Pierpaoli 0], who 
performed a similar investigation using a fiducial two- 
step reionization model instead of the more physically 
motivated fiducial model used here. Their results indi- 
cate that if the fiducial model is not significantly different 
from sudden reionization, assuming sudden reionization 
does not introduce any bias into parameter estimation. 
For a general reionization history, this assumption would 
introduce a bias of 1-3 sigma into values of optical depth, 
amplitude of scalar perturbations and tensor to scalar ra- 
tio. Since we expect the true reionization history to be 
a smooth function of redshift, the conclusion is that for 
near-future CMB data sets, sophisticated modeling and 
reconstruction of the epoch of reionization is unnecessary 
for the accurate determination of other cosmological pa- 
rameters of interest. 
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